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Abstract. 

Thin liquid films on surfaces are part of our everyday life, they serve e.g. as coatings 
or lubricants. The stability of a thin layer is governed by interfacial forces, described 
by the effective interface potential, and has been subject of many studies in the last 
decades. In recent years, the dynamics of thin liquid films came into focus since 
results on the reduction of the glass transition temperature raised new questions on the 
behavior of especially polymeric liquids in confined geometries. The new focus was fired 
by theoretical models that proposed significant implication of the boundary condition 
at the solid/liquid interface on the dynamics of dewetting and the form of a liquid 
front. Our study reflects these recent developments and adds new experimental data 
to corroborate the theoretical models. To probe the solid/liquid boundary condition 
experimentally, different ways are possible, each bearing advantages and disadvantages, 
which will be discussed. Studying liquid flow on a variety of different substrates entails 
a view on the direct implications of the substrate, the experimental focus of this study 
is the variation of the polymer chain length: The results demonstrate that inter-chain 
entanglements and in particular their density close to the interface, originating from 
non-bulk conformations, govern liquid slip of a polymer. 
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1. Introduction 

Confining a material often induces new plienomena tliat cannot be found in bulk. 
Prominent examples are the physics of a quantum dot or the surface reconstructions of 
many crystalline materials. Confinement in liquid flow geometries, however, came into 
focus with the rise of the field of nano- and microfiuidics and with its application in 
lab-on-a-chip devices for transporting, mixing, separating or analyzing minute amounts 
of liquids [U El [3l H] . Recent developments tend to avoid huge external features such as 
pumps to control the flow by designing analogues to capacitors, resistors or diodes that 
are capable to control currents in electronic circuits [5] . Controlling the flow also involves 
a control of the solid/liquid boundary condition. In recent years, various experimental 
methods have been developed in order to characterize the boundary condition to more 
than "no slip" and "full slip" [6l[71[H]. At the same time, theoretical descriptions have 
advanced to include slip effects into the classical thin film equations [H HOI |1T] , which 
now enable a prediction of the influence of shppage to various other situations, e.g. to 
the Rayleigh-Plateau instabihty [121 [131 E] or to spinodal dewetting [15] . Combining 
the experiments with those theories now allows for new insights into conflnement effects 
in liquids, some of which will be reported in this study. 

Liquid slip is characterized by the slip length b and is deflned as the extrapolation 
length of the velocity gradient to zero velocity [16]. Studies of published data of b can 
roughly be divided into two groups, the ones reporting on slip effects in the tens of nm- 
range [TTl dHl [13 [201 [211 [22] and the ones observing large slip effects in the micron or 
the 'plug flow' range [23l[2ll[25l[26l|271|28l[29l[30]. Especially for high molecular weight 
(entangled) polymer melts, large slip lengths were expected according to de Gennes' 
prediction [31] . Moreover, de Gennes proposes the slip length b to increase proportional 
to the third power of the chain length (or the molecular weight, respectively), i.e. 
b = aN^/N^, where a denotes the molecular size and A^^e the entanglement length of the 
polymer. For this study, we especially took care that we will be able to probe the entire 
range of slippage, from no to full slip, and designed the experiments accordingly. 

To probe slip effects in conflned fluids, thin fllms of polymer melts were prepared 
on extremely smooth, well-controlled substrates, typically Si wafers with various 
hydrophobic coatings. To illustrate dewetting. Fig. [T] depicts on the left side a snapshot 
of honey dewetting from a hydrophobic surface and on the right side the lab situation, 
a thin liquid polystyrene fllm dewetting from a hydrophobized Si wafer. For both 
situations, mass conservation is valid and all liquid that once was inside the holes is now 
accumulated in a rim around the holes. As time proceeds, the holes grow. If two holes 
meet, they share a common rim (visible as a straight ridge) that later decays into a row 
of droplets via a Rayleigh-Plateau instability. The stability conditions and the dewetting 
pathways can be characterized according to standard methods [32l |33l IMl [351 [36] . 
Details for the dewetting morphologies, such as nucleation and spinodal dewetting, 
can be found in our recent publications [SS] EZj and references therein. For the 
characterization of the boundary condition in the way it will be presented in the 
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following, the type of dewetting scenario is not relevant; the only precondition is that 
the liquid does not wet the substrate. For all experiments shown in the following, holes 
in the film were initiated by nucleation. 

The aim of the paper is to present the most recent results in a consistent way in 
order to sharpen the conception for slippage phenomena. We therefore concentrate on 
the experimental results and for most of the theoretical details behind the analysis, the 
reader is referred to earlier studies. 



2. Preparation 

As polymers, mostly atactic polystyrene (PS) was used (purchased from PSS, Mainz, 
Germany), but in some cases also block-copolymers are utilized as specified later (c.f. 
section [s]). Atactic PS is preferred as it does not crystaUize, the water uptake is 
negligible and it can be purchased with very low polydispersity and high purity. The 
molecular weight of PS ranged from 5.61 kg/mol (termed 'PS(5.61k)') to 390 kg/mol. 
Polymers are spin-cast onto mica out of a good solvent, which is toluene (Selectipur® 
or LiChrosolv® Line by Merck, Darmstadt, Germany) for PS. The resulting film 
thicknesses are in the range of 100 to 140 nm, unless otherwise noted. Aside from 



single experiments presented in sections 3.2.2 and 3.2.3 PS films exhibiting molecular 



weights M„ > 35.6k were pre-annealed on mica. For the choice of appropriate pre- 
annealing conditions, we made use of experiments presented in Ref. [3H|, where above 
the entanglement length of PS the density of nucleated holes was found to decrease as 
the annealing time increased. Sites of especially high stress might act as nuclei for the 
generation of holes and their relaxation is correlated to the significant reduction of the 




Figure 1. Dewetting films: Honey (left) dewetting from a hydrophobic (OTS-covered 
Si wafer) surface at room temperature and a thin PS(65k) film (right) dewetting from 
a hydrophobized (AF 1600 layer) Si wafer at 140 °C. A typical distribution of nucleated 
holes is observed. 
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hole density. For the largest molecular weight that was used for dewetting experiments, 
i.e. PS(390k), annealing for 3h at 140°C did not lead to a further reduction of the 
hole density, as a plateau value is reached. To safely exclude the presence of internal 
stresses, the aforementioned pre-annealing parameters (3h at 140°C on mica) were 
chosen in a consistent manner prior to dewetting. The thin film was then floated onto 
a MilliQ^^ water surface from where it could be picked up with the sample surface of 
choice. 

Si wafers (Wacker Siltronic AG, Burghausen, Germany and Si-Mat Silicon 
Materials, Kaufering, Germany) have been carefully cleaned using standard procedures 
[5^ and hydrophobized by octadecyl-trichlorosilane (OTS), hexadecyl-trichlorosilane 
(HTS) or dodecyl-trichlorosilane (DTS) self-assembled monolayers (SAMs) using 
standard silanization methods |inillTlll2]- To achieve a hydrophobic coating which is not 
based on silanes, thin Teflon® films were prepared by spin coating a solution of AF 1600 
(Poly[4,5-difiuoro-2,2-bis(trifiuoromethyl)-l,3-dioxide-co-tetrafiuoroethylene], Aldrich) 
in a perfiuorinated solvent (FC- 75®, Acros Organics). Heating the samples above the 
glass transition temperature Tg of the PS film initiates dewetting. 

Tab.[l| gives details of the substrate properties: rms roughness (taken from a 1 /xm^ 
atomic force microscopy scan), water contact angles, static PS contact angle ^y(PS) and 
the surface energies 7sv The latter was calculated from contact angle measurements 
of apolar liquids. Using these parameters together with the dielectric constants and 
refractive indices of all contributing materials, the van-der-Waals (vdW) part of the 
effective interface potential (p can be inferred |13] for each of the hydrophobized 
wafers (see Ref. [37j). Alternatively, experimental access to can be gained recording 
the spinodal wavelength as function of film thickness, as shown in previous studies 
[331 inn [35], [36] . The effective interface potential summarizes all interfacial forces and 
describes the energy needed (or gained) to bring two interfaces to a certain distance. In 
our case, the distance between the two interfaces is the thickness h of the liquid film. 
At values of h where the second derivative of with respect to h, 0"(/i), is negative, 
spinodal dewetting can take place ESI El] • 

Table 1. Properties of AF 1600, OTS, HTS and DTS substrates: thickness d of the 
hydrophobic layer, rms roughness, advancing water contact angle 0adv, water contact 
angle hysteresis A9, surface energy 7sv, static contact angle 9y of PS and absolute 
value of the spreading coefficient S for PS. 
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3. Dewetting as a Probe to Analyze Slippage and Nanorheology 

Fig. [ijright depicts a typical distribution of heterogeneously nucleated holes in a PS (65k) 
film, whereas Fig. |2] shows a temporal series of optical micrographs illustrating the 
growth of a single hole. The impact of the solid/liquid interface to the dewetting 
dynamics is made obvious in Fig. [3^, where the radius i? of a hole is recorded as 
function of time for four different types of hydrophobic layers on top of a Si wafer. For 
the data shown, experimental parameters like PS film thickness (130 nm), dewetting 
temperature (120 °C), molecular weight (13.7 kg/mol), and roughness of the solid 
surface irms < 0.3 nm, c.f. Tab. Q were held constant. Clearly, dewetting on the DTS- 
covered surface proceeds the fastest. The result is surprising, since the driving force 
of dewetting - as characterized by the absolute value of the spreading coefficient S, i.e. 
\S\ = |7iv(cos^Y - 1)1, is largest for AF 1600 and lowest for HTS (c.f. Tab. Q. 
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Figure 2. Growth of a hole in a PS(13.7k) film dewetting from a HTS-covered Si 
wafer at 120 °C. Time interval between subsequent optical images (starting top left 
and ending bottom right): 80s. 
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Figure 3. a) Hole growth dynamics in a thin PS(13.7k) film on various hydrophobic 
substrate surfaces (DTS, HTS, OTS and AF1600) at 120 °C. Here, <=0 is defined 
for each experiment by the first observation of the hole, b) AFM rim profiles (c.f. 
section 3.2.2) of thin PS(65k) films recorded at the same hole radius on DTS, OTS 



and AF1600 substrates. 
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Comparing identically prepared liquid films dewetting from different hydrophobic 
substrates highlights the fact that the fiow properties and the hydrodynamic boundary 
condition at the solid/liquid interface are the key to understand the thin film dynamics. 
This entails the structural and dynamic properties of the interfacial region imposed by 
the substrate and, moreover, also the role of the liquid itself. The latter issue is the 
main subject of this study: We will focus on one hydrophobic substrate, AF1600, and 
systematically vary the liquid properties, in particular the length of the polymer chains 
as given by the molecular weight of the melt. 

Dewetting experiments allow the characterization of the slip length either by 
analyzing the dewetting dynamics [271 EH] or by evaluating the shape of a dewetting 
front [211 |25l |26l |30]. As illustrated in Fig. [3}d, the shape of a liquid rim (resulting 
from identical films and recorded at the same hole radius) is sensitive to the type of 
substrate and, consequently, the boundary condition [29]. Again, different substrates are 
qualitatively compared for the identical liquid (molecular weight, viscosity). A recent 
review describes comprehensively the two mentioned methods to gain the slip length 
and gives detailed views on slip effects in thin polymer films [44j. The main advantage 
of deriving the slip length from the dynamics and the morphology of dewetting films is 
that the method is non-invasive: no external shear has to be applied, no colloidal probe 
or tracer particles interfere with the liquid. 

3.1. Characterizing the Growth of Holes 

From theory [iGj HH HH], one expects a linear growth of the hole radius R with time for 
the no-slip situation and R(t) oc with a = 2/3 in the case of full slippage. Recording 
R{t) curves in holes of dewetting polymer films (c.f. Fig. |2]), a is readily determined by 
fitting an algebraic power law, i.e. R (x {t — to)"- Alternatively, in a double-logarithmic 
plot, the R{t — to) curves are fitted by a line of slope a, c.f. Fig. |4] [55]. The figure 
exemplarily illustrates a set of typical R{t) measurements in various PS films supported 
by AF 1600 substrates. Note that t = marks the time where the hole was first detected 
by optical microscopy; for the evaluation of a, the 'real birth time' of a hole, to, is a 
fit parameter that is determined for each R(t) curve separately. Subsequently, the data 
are evaluated according to the analytical model and analysis technique presented in 
Refs. [271 [28], which is based on the superposition of two dissipation mechanisms [50] : 
viscous friction within the liquid and friction at the solid/liquid interface (slippage). 

3.1.1. Evaluation of the Dewetting Exponent In the following, the molecular weight 
of the polymer melt and the dewetting temperature were systematically varied in 
dewetting experiments. The values for the exponent a given in Fig. |4] already show 
that a does not change significantly with temperature. The effect of the variation 
of M„ illustrates Fig. [5] With increasing M^, a undergoes a transition from high 
values around 0.9 (close to 1) to low values around 0.7 (close to 2/3). The transition 
occurs around the molecular weight of the polymer, where the molecules exceed the 
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critical molecular weight Mc for chain entanglements. As soon as chain entanglements 
become considerable and slippage starts to increase, the dewetting exponent changes 
significantly. This finding is consistent with the theoretical expectation obtained from 
simple energetic considerations, which predict R oc t and R oc t^/^ in the limiting cases of 
no-slip and full-slip, respectively SZl SH] • Deviations from the analytical predictions, 
such as an exponent for the no-slip case being rather close to 0.9 than 1, have also been 
previously reported for numerical simulations [lU] . 

The experimental result, however, visualizes at the same time the limitations of 
this approach (c.f. Fig|5|: Although the slip length varies by orders of magnitude, a 
does not show significant variations except for the transition from weak to strong slip 
(note that the slip length is orders of magnitude larger as compared to the typical film 
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Figure 4. Hole radius R versus time t in a linear (upper left) and logarithmic 
representations obtained from optical micrographs for PS(18k), PS(51.5k) and 
PS(186k) films at different dewetting temperatures. Solid lines represent fit curves 
of the algebraic growth function R oc {t — to)" to the experimental data. Note that a 
shift of the time scale is necessary (the condition R{t = 0) = has to be matched) to 
obtain a linear relation in the double logarithmic representation. 
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thickness). Consequently, drawing further conclusions solely from measurements of the 
dewetting exponent a with regard to slippage and the slip length in particular is not 
appropriate. In anticipation of the next sections, the data for a shown in of Fig. [5^ 
can be related to the slip length b, as depicted in Fig. [5]d. The exponent a is plotted 
as a function of the slip length, which was obtained via the hole growth analysis of the 
R{t) data. The linear superposition of energy dissipation due to viscous flow and due to 
friction at the solid/liquid interface marks the underlying assumption for this method 
[50] and will be introduced in the following. 



3.1.2. Evaluation of the Slip Length by R{t) Experiments Following our description 
as published in Refs. [271 EB], the contact line velocity can be expressed by two 
contributions, resulting from viscous flow and Vs representing liquid slip at the 
solid/liquid interface, given by: 

V = v. + v. = cW-^+^-^- (3.1) 

and the rim width 

w = aV^/R , (3.2) 

where Hq stands for the film thickness and C^{9) denotes the constant of 
proportionality characterizing the fiow field in the vicinity of the three-phase contact 
line (receding contact angle 9). Cg represents the constant of proportionality according 
to conservation of mass and can be determined via in situ AFM measurements (c.f. 
section 3.2.1), while simultaneously recording the evolution of the rim width w and the 



radius R of the hole during dewetting. By plotting the dewetting velocity V (given by 
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Figure 5. a) Dewetting exponent a for a series of dewetting experiments using 
different molecular weights M^. b) Dewetting exponent a versus the slip length b 
as obtained from the hole growth analysis (see section 3.1.2). Colors and symbols of 
data points for different molecular weights M„ in a) correspond to those used in b). 
Dotted lines represent guide to the eyes. 
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the first derivative dR/dt of the hole radius R{t)) versus l/y/R, the viscous and shp 
velocity contributions can be separated. 

Fig. [g] exemplarily illustrates two typical V{l/y/R) curves for PS(65k) films 
exhibiting different film thicknesses on AF 1600 substrates. The film thickness impacts 
the slope of the curves, but does not act on the y-axis intercept Vq. For small hole 
radii, i.e. for large l/\fR, there are data points deviating from the linear regime. This 
deviation can be attributed to the presence of the pre-mature stage of hole growth, 
where the rim of the hole is not fully developed and therefore self-similarity of the rim 
profile, a precondition for the validity of the theoretical description, is not given [271 EI] • 

From the slope of the V{l/\fR) data representation, the slip length h can be 
calculated via Eq. (3.1). Data for the viscosity 77 were extracted from independent 
viscosimeter measurements |52j by extrapolating these data [S3] according to the 
Williams-Landel- Ferry (WLF) equation [Ml 155] . Cg is calculated from AFM 



measurements of the rim width w (see section 3.2) and the film thickness in 
combination with an optical measurement of the hole radius R after quenching the 
sample to room temperature [56] . 

Fig. [7] illustrates the impact of the molecular weight of the polymer melt on the 
slip length h for PS on AF 1600 substrates: Below the critical molecular weight 
for chain entanglements, slip lengths between 10 and 100 nm are found. Above M^, 
the slip length increases with and /3 = 2.9. The values for b{M^), however, are 
systematically lower compared to the expectation by de Gennes (i.e. b = aN^/N^, see 
Ref. [31]) using bulk literature values, namely a = 3 A [17] and Ng = 163 [5^. As will be 




Figure 6. Dewetting velocity plotted versus the reciprocal square root of the 
corresponding hole radius, i.e. V{1/VR). Experiments of PS(65k) dewetting from 
AF 1600 for different film thicknesses at the same dewetting temperature of 130 °C and, 
consequently, identical viscosities. Different slopes of the linear fit curves (solid lines) 
can be found, the y-axis intercept vq, however, is identical within the experimental 
error for all experiments. The dashed lines illustrate the corresponding 95% confidence 
bands of the fit curves. 
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discussed in section 3.2.3 these results are in good agreement with the corresponding 
data for b obtained with the help of the rim profiles of the same holes. 
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Figure 7. Double-logarithmic plot of the slip length b as function of the molecular 
weight A/w of the PS melt; b is obtained from hole growth analysis for single 
experiments (stars) with PS dewetting from AF 1600 substrates. The solid black line 
represents the model by de Gennes [31] using bulk A'^e=163, whereas the dashed line is 
a fit to the experimental data. For reasons of comparison, the solid blue line represents 
the fit to the slip length data (not shown) obtained from the rim profile analysis (see 
section 



3.2.3). 
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Figure 8. Double-logarithmic representation of the y-axis intercept tip, normalized 
by the absolute value of the spreading coefficient versus the reciprocal viscosity 
1/77 (taken from the solid curves in Fig. 14). Mean values of different experiments are 
plotted for different molecular weights below the critical molecular weight M^- 
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3.1.3. Evaluation of the Viscous Velocity Contribution in R{t) Experiments The 
position of the y-axis intercept of the extrapolated hnear fit to the experimental 
V"(1/a/R) data, c.f. Fig. joj is called vo- By plotting fo/|5'| versus l/rj for different M„ < 
Mc, the predicted linear relationship between the viscous velocity contribution Vy = vo 
and the reciprocal viscosity 1/r] is reproduced. This is illustrated in Fig. |8]and represents 
an important consistency check for the applied model. Consequently, the constant of 
proportionality Cy{9) can be calculated from this master curve for unentangled PS 
dewetting from AF 1600 substrates and is determined to be C^iO) = 0.039. Above 
entanglement, slippage changes by orders of magnitude, which has an impact on the 
flow field close to the contact line (c.f. Refs. [IHISZ]). Thus, although still independent 
of the rim size, Cy does depend on the molecular weight for My^ > Mc- 

3.2. Characterizing the Rim Profile 

A rim is theoretically described by a disturbance of a liquid, that decays exponentially 
into the prepared film thickness. Applying thin film models [TOl [26] , the analysis of the 
decay allows for a determination of the slip length, as we have shown in earlier studies 
[211 |25l EHl [29]. The rim profiles were recorded by AFM [57]. The studies shown here 
focus on PS melts of various molecular weights on top of AF 1600. The first part of 
the present section covers the investigation of the temporal and spatial evolution of the 
shape of a liquid front by in situ AFM measurements. The second part deals with 
the evaluation and interpretation of ex situ liquid rim profiles captured at a fixed hole 
radius. The latter experiments and their analysis in view of the slip length have been, 
in part, published in Ref. [3D] and are extended as well as discussed in detail in the 
following sections. 

3.2.1. In situ Measurement of the Evolution of the Rim Shape To probe the temporal 
and spatial growth of the liquid rim, subsequent in situ AFM images during dewetting 
on a high-temperature heating plate were recorded in tapping mode"""^. Cross-sections 
obtained from a temporal series of AFM scans (for a time interval At = 640 s) for 
PS(13.7k) dewetting from AF 1600 at 115 °C are shown in Fig. [9^. The rim clearly 
grows in height and width as dewetting proceeds. Nevertheless, the form of the profiles 
at different times (and therefore different volumes) is very similar. All profiles exhibit a 
substantial trough (see inset of Fig. [9^) on the "wet" side of the rim, where the height 
profile merges into the unperturbed film. 

In order to test the self-similarity of the form of the rim, the rim width w (the 
distance between three-phase contact line and the point, where the rim height has 
dropped to 110% of the initial film thickness) is plotted as function of the maximum 
rim height H, c.f. Fig. |9]d. The linear dependence indicates that self-similarity can safely 
be assumed and w oc \/R (where R denotes the radius of the hole) is obtained from the 
conservation of volume for growing holes (see Ref. [B]). 
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3.2.2. Ex situ Measurement and Analysis of the Rim Shape at Fixed Hole Radii By 
AFM, the rim profile can be recorded at elevated temperature, as described before, 
when the PS is above the glass transition temperature and liquid, or, for simplicity of 
handling the samples, at room temperature {ex situ). The quenching from the liquid to 
the glassy state has in no case lead to differences in the form of the rim profile. It can be 
analyzed as described before [29] (to be more precise, we have used in the following the 
third-order Taylor expanded model based on the full Stokes equations since it applies 
to a non-restricted range of slip lengths, in contrast to e.g. the strong-slip lubrication 
model). For reasons of comparability, all profiles were recorded at a hole radius of 12 fim 
to ensure that approximately the same amount of liquid is collected in each of the rims. 

Rim profiles exhibit two distinct morphological shapes: a) an oscillatory and b) 
a monotonic decay. Exemplarily, Fig. 10 shows typical portions of oscillatory and 
monotonically decaying rim shapes and the corresponding fit functions (exponentially 
damped sinusoidal oscillation and exponential decay, respectively) to the experimental 
data for various molecular weights. The oscillatory decay, which is characterized by 
a significant trough where the rim merges into the unperturbed film ("wet side"), is 
illustrated in Fig. 10 1. In some cases, the experimental data for this part of the rim 
profile even exhibits a slight overshoot as compared to the initial film thickness (see 
inset of Fig. [lO^). With increasing M„, the trough becomes progressively smaller (see 
Fig. [Io|3). Finally, it vanishes and the rim exhibits a monotonically decaying profile (see 
Fig. [Io|;). If two exponential decay lengths can be separated from each other (which 
is always the case for a damped oscillation where the decay lengths are a pair of two 
complex conjugate), the quantification of the slip length b and of the capillary number 
Ca (and also of the viscosity rj by measuring the instantaneous dewetting velocity s 




distance ^ from contact line / nm maximum rim heiglit H / nm 



Figure 9. a) In situ measurement of the temporal evolution (At — 640 s) of the rim 
growth of a 120 nm PS(13.7k) film at 115 °C on AF 1600. All profiles clearly exhibit an 
oscillatory rim shape on the "wet" side of the rim (illustrated by a selection of close-up 
views shown in the inset) . b) Evaluation of the corresponding rim profiles with regard 
to the rim width w and the maximum rim height H. 
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prior to quenching the sample to room temperature) is possible. The determination of 
b in the case of a single exponential decay, however, necessitates the knowledge of Ca 
and, therefore, also of the viscosity t]. The theoretical models and the application of the 
fit functions to the rim profile data are discussed in detail in Ref. [29]. As illustrated 
in Fig. lOi-c and especially Fig. 11 1, profiles show progressively flatter declines for 
increasingly larger molecular weights M^: The rim significantly widens at the expense 
of a decrease in maximum rim height, an observation also reported in Ref. [51]. For a 
given Mw, however, no systematic dependency on the dewetting temperature and, thus, 
on the viscosity of the polymer melt can be recorded (see Fig. 11 d and Fig. 12) [59] . 

To check the influence of possible internal stress in the polymer film due to the 
preparation process [38] [601 [61], the films prepared on mica were annealed prior to the 
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Figure 10. Experimentally obtained AFM images of the "wet side" of rim profiles for 
thin PS films of different molecular weight on AF 1600. a) PS(13.7k) and PS(5.61k) 
(see inset) at 110 °C represent typical rim shapes below the critical molecular weight 
Mc for chain entanglements, b) PS(65k) at 140 °C and PS(81k) at 130 °C (see inset) are 
clearly above M^. c) PS(125k) at 130 °C and PS(271k) at 155 °C (see inset). The solid 
red lines indicate the fits to the experimental data according to a damped oscillation 
or a single exponential decay. 
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transfer to the AF 1600 substrate (c.f. section [2j). A comparison of profiles with and 
without pre-anneahng was exemplarily tested for two molecular weights, PS (65k) and 



PS(125k), and is shown in Fig. 12, For our systems, we could not detect a systematic 



influence of pre-annealing on the form of the rim profile for a fixed molecular weight 



(see also section 3.2.3) 



3.2.3. Evaluation of the Slip Length Fig. [T3| depicts the results of the slip length b 
in view of the molecular weight in a double-logarithmic representation: Below the 
critical molecular weight for chain entanglements M^, the slip length reaches values up 
to 150 nm at maximum. For > Mc, the slip length strongly increases with the 
molecular weight or the length of the polymer chain, respectively. To be more precise, a 
scaling exponent of 2.9(2) between the slip length and the molecular weight, i.e. b oc 
with /3 = 2.9(2), is found. This result is in excellent agreement with the outcome of the 



hole growth dynamics (see section 3.1.2). The values for b are systematically lowered 
compared to the expectation by de Gennes (i.e. b = aN^/N^, see Ref. [31]) based on 
literature values for the molecular size a = 3A [17] and the bulk entanglement length 
Ne = 163 [55]. The pre-factor a/N^ is obtained from the linear fit to the experimental 
data in the double-logarithmic representation: For fixed a = 3A, the rim profile analysis 
leads to = 517, which is a factor of 3.2 larger compared to the bulk entanglement 
length of 163 monomer units [62] . 




Figure 11. a) Examples of rim profiles on AF 1600 for different molecular weights 
Afw Increasing Af^ results in a transition from oscillatory to monotonically decaying 
rims and to successively less steep decays into the unperturbed film, b) Rim profiles 
for PS(18k) as a representative for the typical situation below the critical chain length 
for the formation of entanglements. Dewetting temperature and, thus, viscosity of the 
polymer film have no significant impact on the rim profiles (see also Fig. 12) ,59j . 
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3.2.4- Checks of Consistency and Validity The general consistency of this method to 
extract the shp length and the capillary number from liquid rim profiles was checked for 
PS(13.7k) dewetting from OTS and DTS substrates in our previous publications [261129]: 
The viscosity (which was determined from the capillary number and the dewetting 
velocity) turned out be in excellent agreement with independent (bulk) viscosimeter 
data. In situ measurements did prove that the slip length was independent of the radius 
of the holes (the range between 2.5 and 17 fim was tested). Dewetting experiments for 
various film thicknesses (between 50 and 230 nm) also lead to consistent results. The 
method, consequently, yields reliable results for the slip length and the viscosity of 
dewetting thin polymer films on the silanized Si wafers. In the following, the identical 
consistency and validity checks are performed for the dewetting experiments on the 
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Figure 12. Testing the role of pre-annealing (films with or without pre-annealing for 
3h at 140° C on mica prior to the transfer to the AF 1600 substrate) and dewetting 
temperature for the shape of the rim profile on AF 1600 for PS(65k), PS(125k), 
PS(186k) and PS(390k) [Sg. 
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Determination of the PS Film Viscosity For AF 1600 substrates, the capillary 
number Ca is extracted from the analysis of the experimentally recorded rim profiles 
up to PS (81k). The knowledge of the surface tension of the polymeric liquid 
(7iv(PS) = 30.8 mN/m, [63]) and the experimental determination of the actual dewetting 
velocity s (from optical hole growth measurements) enable the calculation of the 
viscosity rj of the polymer melt via Ca = rjs/'jh for a given dewetting temperature. 
Dewetting experiments were performed for various molecular weights of PS melts below 
and above the critical molecular weight Mc for entanglements. The results are given in 
Fig. 14 and Tab. [2j Within the experimental error, good agreement with independent 
viscosimeter measurements and their extrapolation according to the WLF equation 
[52| l53| [SH |55] is achieved. From monotonically decaying rims only one decay length 
can be obtained. Therefore, the capillary number and, thus, the viscosity cannot be 
extracted for PS(lOlk) and larger on the AF 1600 substrate. 



Preparation Conditions and N on- Newtonian Effects In the context of the results of 
this dewetting study concerning the flow dynamics of thin polymer films, a discussion 
of the influence of potential viscoelastic effects is necessary. A source of viscoelastic 
effects are residual stresses caused by the fast solvent evaporation during the spin- 
coating process. To safely exclude this issue, dewetting experiments were typically 
performed after a pre-anneaiing step on the mica substrate (c.f. section [2]). However, 
to test the influence of this procedure, additional experiments were carried out without 
annealing prior to dewetting. Rim profiles of these as-cast films were analyzed and 
compared to pre-annealed samples: No systematic influence of the pre-anneaiing step 
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Figure 13. Double-logarithmic plot of the slip length 6 as a function of molecular 
weight My^. b is obtained from the analysis of rim profiles for single experiments 
(symbols) with PS dewetting from AF 1600 substrates. The solid black line represents 
the model by de Gennes [3T] using bulk 7Ve=163, whereas the solid red line is a fit to 
the experimental data. Additionally, film thickness and pre-anneaiing conditions were 
varied. 
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Table 2. PS on AF 1600: Quantitative comparison of viscosities ?7iim determined from 
capillary numbers which were obtained by the analysis of rim profiles to viscosities 
VWLF calculated from the WLF equation. 
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is found for the typical molecular weights of PS (65k) and PS (125k) used in this study 
(see Fig. 12). The slip lengths determined from these experiments are still roughly 
one order of magnitude lower than expected from literature values for A^^e (see Fig. 13). 
Consequently, for "mature" holes exhibiting fully developed rims, the flow dynamics and 
morphologies are not significantly affected by polymer relaxation processes that could 
potentially originate from thin film preparation and residual stresses. 

Further the question arises whether the flow itself is Newtonian, i.e. whether 
the viscosity is independent of the shear rate within the range of shear rates present. 




100 110 120 130 140 
temperature T / °C 



Figure 14. Semi-logarithmic plot of viscosity 77 versus temperature T for PS melts of 
various molecular weights ilf^ < -^j'^c- The viscosity data is obtained from rim profile 
analysis of PS(5.61k), PS(13.7k) and PS(18k) films on AF 1600. The solid curves are 
based on the extrapolation of a fit to experimental viscosimeter data for PS(12.5k) 
according to the WLF equation. 
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Weissenberg numbers Wi = t'j smaller than 1 assure that viscoelastic effects of the 
polymer melt can safely be neglected [61]. The present experimental configuration, 
given by the dewetting temperatures and the molecular weights of the polymer, results 
in characteristic relaxation times, shear rates, and Weissenberg numbers for which 
estimates can be calculated (see Tab. |3|: Relaxation times r are determined via the 
viscosity rj, obtained from the WLF equation (c.f. Refs. [52|, |53l ISH [55]), and the shear 
modulus G of the polymer melt. Shear rates 7 were estimated as described in Ref. jSU] 
resulting in an upper limit for the corresponding Weissenberg numbers Wi=T'y. As 
indicated in Tab. |3| Weissenberg numbers calculated from experiments on the AF 1600 
substrate do not exceed a value of 0.1. Consequently, the occurrence of non- Newtonian 
effects is not expected for the examined molecular weights and the low shear rates 
(see Tab. |3| in such dewetting experiments. The analysis of the dewetting experiments 
presented in this study can therefore safely be based on the assumption of Newtonian 
flow. 

Impact of the Film Thickness Finally, the question arises whether slippage of a thin 
liquid film, representing fluid motion in a confined geometry, depends on the level of 
confinement. To test the effect of the film thickness on the slip length, an additional 
dewetting experiment (without pre-annealing) was performed using PS(271k) at 140 °C 
with a film thickness of 260 nm. The results are presented within Fig. [13] as the red data 
points: The analysis of the rim profile of this hole gave a slip length of 12.4 /im. This 
is, within the experimental error, in very good agreement with the corresponding slip 

Table 3. Calculated relaxation times t — -q/G and Weissenberg numbers Wi for 
diverse temperatures T and molecular weights of the PS melt dewetting from the 
AF 1600 substrate. The viscosity 77 was obtained from the WLF equation; G ~ 0.2 MPa 
was taken from literature |55j . The maximum shear rates 7niax and the maximum 
values for corresponding Weissenberg numbers Wi were calculated as described in 
Ref. [30]. 
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length of 11.8 /im for a 115 nm thick PS(271k) film. A further hole growth experiment in 
a 230 nm PS(65k) film (after a typical pre-annealing of the film on mica) was recorded 
at 130 °C. The analysis of the corresponding rim profile exhibited a slip length of 73 nm. 
Both experiments using thicker films do not show a significant increase in the slip length. 



4. Discussion 



In literature, slip lengths determined from molecular dynamics (MD) simulations in 
non- wetting conditions are reported to be in the range of 15 molecular diameters [65] . 
Approximating the molecular diameter by twice the radius of gyration Rg of a polymer 
coil (for PS(13.7k): Rg = 3.2 nm) results in values for the slip length of up to 100 nm. 
From this point of view, the experimental results for the slip length on AF 1600 below 
Mc are in good accordance with the theoretical expectation. 

For molecular weights above Mc, b oc is found in accordance with de Gennes' 
prediction. However, using the bulk value for the entanglement length, = 163, slip 
lengths are expected to be roughly one order of magnitude larger than what is found in 
our experiments. The deviations can be explained by a reduced entanglement density 
compared to bulk values. This is quantified by an increased entanglement length of 

iVe=517. 

According to the fact that slip lengths are reproduced for markedly thicker films, 
we conclude that the observed reduction of the entanglement density is not caused 
by the confinement of the film but is localized at the solid/liquid interface. The 
lack of a significant temperature dependence highlights the fact that slippage and 
fiow dynamics are directly related to inter-chain entanglements rather than viscosity: 
The difference in chain length between e.g. PS(51.5k) at T=125°C and PS(271k) at 
T=155°C is responsible for a increase in slip of two orders of magnitude, while the 
viscosity (2 x 10^ Pa s) is identical for the two examples [30]. 



4-1- Reduced Interfacial Entanglement Density and its Implications 

The analysis of the slip length as a function of the molecular weight of the polymer has 
shown that close to the interface, deviations from bulk properties are to be expected. 
Polymer chain entanglements and hence the entanglement length of the polymer are 
related to the packing of the polymer chains in the bulk. Close to an interface, however, 
the chains pack differently. The reduced entanglement density in the vicinity of the 
solid/liquid interface results from a reduction in the volume pervaded by a polymer 
chain at an interface as described by Brown and Russell [BB]: The polymer coil cannot 
cross the boundary at the interface that rather acts as a reflecting plane, as previously 
proposed by Silberberg [HZl EH] and is known as Silberberg's principle. Chain segments 
fold back on the same chain and, thus, the degree of coil-coil interpenetration decreases 



c.f. Fig. 15 ). The reflection process causes a decrease of the radius of gyration normal to 



the plane, but does not affect the radius of gyration parallel to the plane. The polymer 
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coil can be imagined having an ellipsoidal rather than a spherical shape j69]. Close to 
the interface or the wall, the volume pervaded by a given chain length is smaller than 
in the bulk. Brown and Russell [66] suggest that the vertical distance affected by this 
phenomenon can be approximated at most by the radius of gyration Rg of the polymer. 
Hence, the normal bulk entanglement density is expected to be recovered at a distance 
of Rg above the interface. They analytically calculate a reduction of the entanglement 
density or, analogous to that, an increase of the entanglement length by approximately 
a factor of 4, which is in accordance with the experimental results presented here. 

The impact of the density of inter-chain entanglements on the flow profile and 
the slip velocity at the interface of an entangled polymer melt and a solid substrate is 
difficult to resolve on the molecular level. For a simple liquid, true and apparent slip 
can be distinguished according to the velocity of the first as compared to subsequent 
layers of molecules (see e.g. Ref. |44j)- In contrast to this, the motion of polymer 
chains in entangled melts is much more complex, as e.g. De Gennes' reptation model 
shows. To interpret our results, we propose an apparent slip velocity that depends 
on the entanglement density within this interfacial layer. This picture does not take 
into account relative motion on the segmental level within this region and, therefore, 
cannot answer the question whether the monomers or even the chain segments that are 
in contact with the substrate exhibit a non-zero velocity. The apparent slip velocity, 
representing the entire interfacial region, shifts the fiow profile within the liquid film 




Figure 15. Polymer conformation in a bulk melt (upper figure) and at the 
substrate/melt interface (lower figure): Two neighboring coils interpenetrate each 
other and may form entanglements. At the interface, the radius of gyration parallel 
to the substrate Rg(\\) is comparable to bulk i?g, while it is reduced in perpendicular 
direction (i?g(j^)). Consequently, the coil-coil interpenetration decreases resulting in a 
loss of entanglement at the substrate/melt interface. Inspired by Ref. [70] . 
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towards larger velocities. In the strong-slip situation that we are reporting here for 
entangled polymers, the flow profile is linear and its gradient is unaffected by slippage 
as it depends on the viscosity of the (bulk-like) residual film. Fig. 16 illustrates the 
difference between apparent slip for a reduced interfacial entanglement density and for 
a bulk-like entanglement network. 

Recent experimental studies probing the effect of confined polymer chains in thin 
films have also provided strong evidence of deviations from polymer bulk properties 
and a reduced entanglement density: Si et al. report a decrease of the entanglement 
density in thin freestanding films compared to the bulk polymer [71] . They conclude that 
the total entanglement density is constant, but due to the reduced chain overlap near 
interfaces, the proportion of self-entanglements increases at the expense of inter-chain 
entanglements. Rowland et al. recently studied the squeeze flow of entangled polymers 
in confined conditions of films thinner than the size of the bulk macromolecule p2] . Their 
experiments reveal a reduced interpenetration of neighboring chains in favor of self- 
penetration leading to a weakening of the polymer network and a decreased resistance 
to melt flow due to the presence of less inter-chain entanglements. Molecular dynamics 
simulations and corresponding primitive path analysis for polymers studying chain 
conformations and local dynamics corroborate the occurrence of a reduced entanglement 
density near a wall [721 El] • Further studies address the film preparation process as a 
source of non-equilibrium conformations: Barbero and Steiner studied the mobility at 
the surface, i.e. at the liquid/air interface of thin polymer films by an electric field 
technique [75]. They attributed the reduced effective viscosity in the case of as-cast 
films compared to annealed ones to a reduced, non-equilibrium entanglement density 




Figure 16. AFM height profile of a dewetting rim (dewetting velocity V) and (non- 
scaled) illustration of the flow profile (velocity field u) for entangled polymer films. 
The apparent slip velocity (blue arrow), representing the interfacial region at the 
solid/liquid interface as described in the text, is determined by the local density of 
inter-chain entanglements. A reduced interfacial entanglement density, compared to 
the bulk entanglement network, is responsible for a significantly smaller slip velocity 
and a shift of the entire flow profile towards smaller velocities. The slip length b is 
obtained as the extrapolation length of the flow profile to zero velocity. 
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resulting from preparation. 

The effect of an altered entanglement network and, in particular, a reduced 
entanglement density at an interface also impacts polymer/polymer (and their 
interdiffusion) or polymer/air interfaces, for which the hypothesis of the reflection of 
polymer segments is applicable, too. The unique feature of the experimental results 
presented here consists in the fact that sliding friction is sensitive to the proximate 
properties of the polymer at the solid/liquid interface, though the residual bulk film 
is not necessarily affected by confinement. The above mentioned experimental studies 
[7T| [72] relate their findings to the confinement given by the film thickness (and its 
ratio to the spatial dimensions of the polymer chains therein). Here, however, the direct 
and unique impact of the solid/liquid interface upon chain conformations and rheological 
properties of the polymers in the vicinity of a substrate or wall is revealed. The thickness 
of the layer affected by the reduced entanglement density at the solid/liquid interface 
is at most in the range of the radius of gyration [66]. i?g depends on the chain 
length of the polymer and varies for typical PS molecular weights > Mc used 
in this study from 5.2 nm for PS(35.6k) up to 17.1 nm for PS(390k) [76]. Moreover, 
at the polymer/air interface, a reduced entanglement density might be present that - 
potentially - propagates the same distance away from the interface into the polymer 
film. For films thicker than 100 nm, the "bulk region" of the film is even for the largest 
M„ of our study still more than two times larger than the superposition of both effects 
(layers). In some studies, however, Rg is chosen to be larger than the film thickness, 
leading to effects of "overlapping interfacial layers". Hence, these special cases cannot 
be captured via the methods presented here. These cases are especially dependent on 
the preparation procedure and show massive effects upon pre-annealing [SOI • 

The decrease of the entanglement density at the solid/liquid interface is, as 
described before, inferred from the determination of the entanglement length = 517. 
This value corresponds to a molecular weight of the polymer strand between two 
entanglements of Me = 53.8kg/mol. Thus, being significantly larger than the critical 
molecular weight found for bulk situations, this consideration suggests that the 
formation of entanglements at the interface starts at a higher molecular weight as 
compared to the upper residual (bulk-like) polymer film. 

4-2. Quantitative Comparison of Slip Lengths: Rim Profile Analysis vs. Hole Growth 
Dynamics 

Slip lengths obtained from hole growth dynamics agree qualitatively with those gained 
from the analysis of rim profiles: Slip is independent of the molecular weight for 
unentangled films, whereas the slip length increases oc for molecular weights above 

While results are also in good quantitative agreement for low molecular weights, 
however, a significant deviation is found for M„ > M^. data extracted from hole 
growth are typically larger by a factor of about 2-3 compared to data obtained from 
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the rim shapes (c.f. Fig. [7|). This difference is potentially caused by an inevitable 
assumption included in the energy dissipation model of sliding friction: Regarding 
the two dimensional representation of the rim, the lateral extent of slippage in the 
x-direction, i.e. in the direction of dewetting, is approximated by the rim width w. The 
latter is again connected to the slip length: Increasingly larger slip lengths are correlated 
to wider rims. The rim width, however, has been defined for practical reasons as the 
distance between the three-phase contact line and the x-coordinate at which the rim 
height has dropped to 1.1 times the initial film thickness h^. If, however, the "true" 
lateral extent of slippage is significantly smaller than the rim width w, the frictional 
energy is dissipated over a narrower extent and the slip lengths obtained from the model 
for the hole growth analysis (assuming a superposition of viscous flow and slippage, 
c.f. Eq. (3.1)) would be systematically too large. This hypothesis is supported by the 
dewetting experiments on the AF 1600 substrate below the critical chain length for 
entanglements: In case of these solely weakly slipping PS films, no significant difference 
between slip lengths obtained from rim profiles and hole growth dynamics is found. 
Consequently, the significant shift between the results of the two methods is either 
correlated to the change of the rim profile, in combination with the definition of the rim 
width, or, and this is the second possibility to be taken into consideration, to the change 
of the velocity profile according to the presence of slippage. Gaining information about 
the extent of slippage is difficult. In section |5} an experimental system will be introduced 
that might help to shed light on this issue: thin block-copolymer films dewetting from 
hydrophobic substrates. 

Besides the aforementioned methodical issue of the differences concerning the 
determination of slip lengths from the hole growth dynamics and the analysis of rim 
profiles, a further explanation can be hypothesized: Rim profiles are sensitive to the 
local hydrodynamic boundary condition in the region of the solid/liquid interface where 
the rim decays into the unperturbed film. The hole growth analysis, however, is sensitive 
to the entire extent of the rim width w. 

To sum up the above discussion, the systematic differences in slip length originating 
from the two different techniques may be linked to the lateral extent of slippage 
underneath the liquid rim. 



5. Pushing to the Limits - More Sophisticated Dewetting Experiments 

In the previous section, we addressed the question to what lateral extent the area 
underneath the rim is exposed to interfacial flow and in particular interfacial frictional 
forces: Which part of the liquid rim is actually capable to slide over the substrate's 
surface? To explore this issue, we modifled the experimental system and switched to 
liquids with an internal structure, block-copolymers. The idea is to prepare a well-known 
di- or triblock-copolymer microphase-separated thin fllm and then to induce flow by 
dewetting from a hydrophobic substrate. Similar to homopolymer fllms, holes break up 
and grow in time while forming a liquid rim. At distances from the three-phase contact 
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line where the interfacial microphase structure is exposed to shear forces, a structural 
orientation of the cylindrical micro-domains is expected to be induced; the phenomena 
of shear-induced orientation has already been reported by other studies before [771 EH] • 
After quenching the sample to room temperature, the film is detached from 
the substrate and the rim is imaged from underneath with the AFM. Differences 
in the orientational structure of the two phases become visible for different lateral 
displacements from the three-phase contact line [TO]: Fig. [Tt] depicts an AFM (tapping 
mode"'"'^) phase image of a PS/PEP block-copolymer (Kraton® G1701E, volume 
percentage of PS: 28%, Kraton Polymers, Houston, TX, USA) film dewetting from 
a hydrophobized (OTS-covered) Si wafer. A comparison of the width of the oriented 
area, as e.g. determined from Fourier transforms of consecutive small areas, to the 
topographical rim width can now be accomplished. Roughly, the width of the oriented 
area is 1-2 /im. The range corresponds to the width of the rim. (Details will be published 
elsewhere [79].) Ongoing studies focus on the systematic modification of slip properties 
by a variation of the type of substrate, which has shown to be strongly correlated to slip 
of PS homopolymer melts (c.f. Fig. [s]), and further parameters such as the dewetting 
temperature [TH]. In case of a strong-slip boundary condition, the structure close to 
the solid/liquid interface should experience friction rather than significant shear: it is 
rubbed over the surface. This should affect the microphase structure in a different way 
as compared to shearing, which is related to a no- or weak-slip boundary condition. 




Figure 17. Microphase structure (AFM phase image) of a PS/PEP block-copolymer 
film at the substrate/polymer interface as recorded from underneath. The image shows 
the vicinity of an hole (bright part on the left side), which resulted from dewetting of 
the film from a hydrophobized (OTS) Si wafer at 140°C. 



Slippage and Nanorheology of Thin Liquid Polymer Films 



25 



6. Conclusion and Outlook 

The challenge of explaining facets such as sliding friction and the fluid dynamics of 
polymers at the solid/liquid interface on the molecular level is a main task in the field 
of micro- and nanofiuidics. The aim of this study was the characterization of the flow 
dynamics of thin polymer films in view of the hydrodynamic boundary condition at 
the solid/liquid interface. Dewetting experiments are a powerful tool to study the 
mechanisms of friction and the rheological properties of polymers on the nanoscale. 

Hole growth dynamics and rim profile analysis allow for a quantification of the 
impact of the molecular weight and, in particular, of chain entanglements on the sliding 
properties of polymers: On smooth, hydrophobic AF 1600 surfaces, a weak-slip situation 
is found for PS melts with molecular weights below the critical molecular weight 
for entanglement formation. This result is independent of the molecular weight in this 
regime. However, slip is strongly enhanced above M^: The slip length scales with the 
third power of the molecular weight of the polymer, as predicted by De Gennes. At the 
interface itself, the density of chain entanglements is found to be reduced by a factor of 
3 to 4 compared to the bulk, caused by the different packing of polymer chains at an 
interface according to Silberberg's principle. 

The results are of extraordinary importance in view of a series of applications such 
as e.g. polymer extruders or, in general, microfiuidic devices. As pointed out by Bocquet 
and Barrat |8], the mean velocity of a pressure drop flow in a cylindrical channel is 
increased by a factor 1 + 86//i as compared to the classical no-slip boundary condition, 
where h is the slip length and h the thickness of the channel. Thus, the permeability 
of channels and also of porous media is highly increased by slippage. An increase in 
the slip length for longer polymer chains by two orders of magnitude, as presented in 
this study, would also impact the mean velocity in sub-micron channels by orders of 
magnitude. 

In addition to the influence of polymer chain length on slippage discussed in 
this paper, we shortly reported on significant dynamical and morphological differences 
purely resulting from the choice of different hydrophobic substrates (c.f. Fig. [s]). In 
terms of drawing a general molecular picture, not only the degree of entanglement 
but also structural and dynamic properties on the monomer and segmental level 
imposed by the substrate have to be discussed. As presented in Tab. [T], characteristic 
substrate properties such as roughness (based on AFM measurements), macroscopically 
determined wetting properties and surface energies, however, reveal no clear correlation 
of one of these parameters to slippage. For a discussion of the interfacial structure on 
the molecular level and its implications on the dynamical properties of the liquid, the 
reader is referred to our recent study ||80j : Neutron and x-ray reflectometry are powerful 
experimental tools to explore the molecular structure of the solid/liquid interface and 
add new insights explaining the differences in the sliding properties of PS on different 
substrates. Currently, also dewetting experiments are extended: The aforementioned 
methods and theoretical models can also be applied to the spatial and temporal evolution 
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of thin poly(methyl methacrylate) (PMMA) films. The results might shed further 
light on the link between the monomeric structure of the polymer and slippage at the 
solid/liquid interface. 

Aside from the analysis of the hole growth dynamics and the morphology of the rim 
surrounding the hole, recent experiments and theoretical studies also focus on the liquid 
rim instability as an indicator for the slip/no-slip boundary condition [121 1131 [IH l8T] . 
This instability is comparable to the Rayleigh-Plateau instability emerging on a surface 
of a cylindric jet of water flowing off a tap. The resulting surface corrugations can be 
optically observed and the jet finally decays into discrete droplets. The morphology 
of the instability of retracting straight fronts in thin liquid films is governed by the 
hydrodynamic boundary condition: While under no-slip conditions symmetric liquid 
bulges (co-moving with the front) are found, slippage causes the formation of asymmetric 
finger-like structures. This phenomenon can be optically observed and represents a 
further way to access slippage in thin polymer films. Further characteristics of the 
instability such as the evolution of the wavelength of the instability along the receding 
front (on lateral scale and perpendicular to the surface plane) are subject of ongoing 
investigations [ST] . 

We hope that our results will motivate further theoretical studies and simulations, 
as this phenomenon could be of high relevance for other experiments such as e.g. colloidal 
particles moving through a liquid medium. 
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